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Results from the most extensive study of the time-evolving dust structure 
around the prototype "Pinwheel" nebula WR 104 are presented. Encompassing 
11 epochs in three near- infrared filter bandpasses, a homogeneous imaging data 
set spanning more than 6 years (or 10 orbits) is presented. Data were obtained 
from the highly successful Keck Aperture Masking Experiment, which can recover 
high fidelity images at extremely high angular resolutions, revealing the geome- 
try of the plume with unprecedented precision. Inferred properties for the (un- 
resolved) underlying binary and wind system are orbital period 241.5±0.5 days 
and angular outflow velocity of 0.28±0. 02 mas/day. An optically thin cavity of 
angular size 13.3 ± 1.4 mas was found to lie between the central binary and the 
onset of the spiral dust plume. Rotational motion of the wind system induced 
by the binary orbit is found to have important ramifications: entanglement of 
the winds results in strong shock activity far downstream from the nose of the 
bowshock. The far greater fraction of the winds participating in the collision 
may play a key role in gas compression and the nucleation of dust at large radii 
from the central binary and shock stagnation point. Investigation of the effects 
of radiative braking pointed towards significant modifications of the simple hy- 
drostatic colliding wind geometry, extending the relevance of this phenomena to 
wider binary systems than previously considered. Limits placed on the maxi- 
mum allowed orbital eccentricity of e ^ 0.06 argue strongly for a prehistory of 
tidal circularization in this system. Finally we discuss the implications of Earth's 
polar (z<16°) vantage point onto a system likely to host supernova explosions at 
future epochs. 

Subject headings: stars: individual(WR 104), stars: mass loss, instrumentation: 
interferometers, stars: imaging, stars:winds, outflows 



Introduction 



The first well-resolved images of the dusty, IR-luminous Wolf-Rayet star WR 104 re- 
vealed an elegant plume stretching hundreds of AU from the bright co re and following a 
trajectory closely matched to an Archimedeian spiral (jTuthill et al.lll999l ). The origin of this 
geometry is simple and highly intuitive: material is embedded within a uniformly-expanding 
spherical wind driven from the hot stars at the heart of the system, yet the orbital motion 
of a central binary causes a rotating wake embedded within the flow. These images gave 
immediate confirmation of the colliding-wind binary nature of dust production in WR 104, 
and at the same time presented astronomy with a new type of object for study: a "Pinwheel" 
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nebula. 



Although we found no specific predictions in the literature for such a spiral dust plume 
in dusty Wolf-Ra yets such as WR 104, the same combination of orbital and outflow motions 
led iKuiperl (1194 ll ) to predict exactly this geometry for the circumstellar environment around 
the eclipsing binary (5 Lyr. In fact use of th e term "Pinwheel" to describe such a structure 
can be found as early as 1950 (IStruvdll950l ). and theor etical study of colliding-wind binary 
syste ms has since developed an extensive literature (e.g. lStevens et al.lll992uWalder fc Folini 
2000l . and refs therein). The basic picture of colliding- winds mediating dust production was 
formulated to explain the episodic dust-producer WR 140, which be came something of a 



'Rosetta S tone' system in important papers establishing the model (jWilliams et al.l 11990 



Usovl Il99ll ). It is also interesting, although perhaps tangential to the work of this paper, 
to note that Archimedian spirals may arise in quite distinct mass-loss e nvironments such as 
the b inary reflex-motion induced structure recently seen in AFGL 3068 OMauron fc Huggins 
2006h . 



Observational confirmation of the pinwheel nebula around WR 104 has not yet been 
attained by adaptive optics imaging due to the requirement for a stable, well characterized 



point -spread function enabling recovery of structure at the diffraction limit (IRajagopal et al. 



20041 ). However, a number of other techniques have been used to confirm the presence of the 
spira l plume including high spatial resolution lunar occultations (IMondal fc Chandrasekhar 
20021 ) . variations in t he visible lightcur ve of 241 days in accord with the orbital period derived 
from interferometry ( Kato et al. 2002h . together with exten s ive follow-on observa tions with 
aperture masking interferometry (ITuthill et al. I2002L 120031 : iMonnier et~aD l2007h . Reports 
of binarity in the WR 104 system (jCrowtherl 119971 ) appeared shortly befo re the original 



discov ery of the spiral, with the spectral classification for the object given in Ivan der Hucht 



(120011 ) as WC9d + B0.5V (+VB), although given the uncertainties we refer to the second 
component spectrum (actually brighter by a factor of ~2) as OB. Detailed radiative transfer 
modelling, confirming the appearance and basic observational properties of the pinwheel, 
has been performed by lHarries et al.l (120041 ) . 



The discovery of similar structures in an expanding list of host systems has served to 
motivate further inte rest in these objects. The second such system to be found, WR 98a 
( IMonnier et al.lll999l ). showed a similar morphology implying approximately matched basic 
properties (in terms of orbital period and windspeed) to the prototy pe WR 104. However, a 
more diverse population is also coming to light including WR 112 (IMarchenko et al.ll2002l ) 
which appears to have a significantly longer 25 yr binary period, and the high orbital eccen- 
tricity 8 yr system WR 140 whos e periodic periastron passages have been shown t o throw off 



a remarkable series of dusty arcs (IMonnier et al.ll2002 



Marchenko fc Moff atll2006l ) . Recently, 
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a dramatic association of at least two and probably five such spiral sys tems was discovere d 
near the Galactic Center at the heart of the massive Quintuplet Cluster ( ITuthill et al.ll2006l ). 

The status of Wolf-Raye t s as candidate progenitors to type Ib/c supe r novae (reviews 



can be found in Bet he 1990; Vanbeveren et al. 



19981 : IWooslev et all 120021 ) has generated 



considerable interest in their mass loss histories. In particular, structures within the cir- 
cumstellar environment, into which the supernova blast wave expands, have the capacity 
to interact w ith the ejecta mod ifying the lightcurve and observed properties of the explo- 
sion. Indeed, iRyder et al.l (120041 ) observed a residual fluctuating or rippled signal, obtained 
when subtracting a spherical best-fit model from the observed radio lightcurve of super- 
nova 200 lig. The most likely interpretation for these modulations, they argue, arises as the 
supernova shock encounters periodic density enhancements embedded in the circumstellar 
environment by a colliding-wind binary pinw heel system. With the ra pidly solidifying links 
between supernovae and gamma-ray bursts (IWooslev fc Blooml 120061 ) . the mechanisms by 
which such mass-loss signatures may become encode d within the lightcurves of the most 
energetic explosions in the cosmos are being explored (IMoran fc Reichartl 120051 ). 

In this paper, we present the results of the most extensive diffraction-limited imaging 
study yet performed on the prototype pinwheel system, WR 104. Observations encompass 
three different near-infrared filter bandpasses and eleven distinct epochs covering more than 
six years (almost ten orbital periods) of the system. With the structure of the time-evolving 
plume revealed in unprecedented detail, we are able to address key outstanding questions re- 
garding the way in which the dust formation is mediated by the colliding winds of the binary. 
Multi-wavelength observations are helpful in teasing apart the complex three-dimensional ge- 
ometrical, thermal and illumination effects which go into creating the phenomena we see as 
a pinwheel nebula. 



2. Observations 



Monitoring of a number of dust-forming Wolf-Rayet stars was performed as part of our 
aperture-masking interferometry program begun at the Keck 1 telescope in 1996. This has 
produced astrophysical results spanning a r ange of different are a s within contempor ary stellar 



physics including studies of young stars (IDanchi et al 



pulsating stars and giants (ITuthill et al.l 12000 



2001 



Monnier et al 



shrou ds in proto-planetary nebulae and transition objects ( IMonnier et al.ll2000t ITuthill et al. 



Tuthill et all 120021 ): evolved 
2004T); and dusty mass-loss 



20021 ). A full description of the experiment including a discussion of the concept ual principles 



and si gnal-to-noise considerations underlying masking interferometry is given in ITuthill et al. 



( I2000l ) while further discussion of systematics and seeing induced errors can be found in 



- 5 - 



Monnier et all (12004 ). 



Observa tions presented here all u sed a partially-redundant annular mask geometry first 
suggested by lHaniff fc Buscherl (119921 ) which was found to deliver robust imaging for targets 
as faint as nix ~ 6 mag with complete Fourier coverage out to the 8 m maximum baseline 
passed by the mask. The detector used for this work was the NIRC camera with image 
magni fying front-optics mod ifying the plate scale so as to sample the diffraction-limited 
beam jMatthews et alJll996h . A relatively rapid-exposure (T int = 0.14 sec) and high data 
volume (data cubes containing 100 frames) speckle observing mode was utilized, although 
this was not as well developed as it could have been and suffered from an inherent low duty 
cycle of at best 20% set by the camera electronics. 

Observations of science targets were interleaved with those of nearby unresolved point- 
spread function (PSF) reference stars, and the extraction of calibrated visibility a nd closure- 
phase data followed established practices within speckle or masking interferometry (jTuthill et al 
20001 ) . All V 2 and closure phase data product s were stored in t he FITS-based Optical In- 
terferometry data exchange format (OI-FITS iPauls et al.ll2005l ). and are available on re- 
quest from the authors. Images were r ecovered using a number of different numerical meth- 
ods including the CL EAN algorithm (H ogboml Il974f) and with severa l maximum entropy 
JGuII fc Skirling! Il984h based methods Tsivial Il987i llreland et al.l l2006h . With some varia- 
tion in levels of fidelity and resolution obtained, all methods produced similar structures to 
those presented in the f ollowing se ctions, and hereafter we show exclusively the results of 
the VLBMEM package <jSivialll987h . 



A journal of masking observations of WR 104 is given in Table [IJ while information about 
the three interference filters within the infrared H, K and L bands utilized in the study is 
given in Table [2J A number of different PSF reference stars were observed, all of which are 
relatively nearby (typically within about 5° on the sky) and tailored to give roughly similar 
received counts under various observing configurations. The primary calibrator stars utilized 
were SAO 186681, SAO 186841, 14 Sgr and HD 165813 although occasionally additional stars 
were used for confirmatory checks. 

In addition to their scientific merit, it is worth noting that the images presented in 
later sections are a strong testament to the power of aperture masking interferometry in the 
recovery of high fidelity images. In particular, the pinwheel around WR 104 resolved here 
possesses strongly asymmetric structure at intermediate dynamic range to the bright core. 
Deconvolution to obtain at least the full diffraction limit is essential (and ideally, some super- 
resolution beyond it), as the structure recovered is formally only a few resolution elements 
across. Seeing conditions over the twelve observing nights found in Table [1] varied over a 
wide range from extremely poor to very good. Robust, high-fidelity images were recovered 
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for all epochs and in all niters observed. 

This is in contrast to the failure to recover significant reliable structure from the on e 
attempt (we are aware of) to image this target with adaptive optics (IRajagopal et al.ll2004j ). 
The presence of asymmetric and time-varying extended structure at high resolutions which 
is nevertheless completely predictable in morphology makes WR 104 almost ideally suited 
as a "test pattern" target against which imaging experiments may be benchmarked, and we 
look forward to confirmatory results from new generations of AO instruments. 



2.1. K-band images 

Images recovered from data taken at 11 separate epochs in the CH4 filter are presented in 
Figured! Although there are 12 observing dates given in Table (H the first two of 1998 Apr 14 
and Apr 15, being separated by only one day, did not give significantly different images and 
these data have been averaged together for the remainder of this paper. Experience has 
taught that images recovered from filters within the infrared K-band usually deliver the 
highest fidelity. The PAHCS filter at longer wavelengths has significantly lower angular 
resolution with a corresponding loss of detail. On the other hand, the H filter at shorter 
wavelengths suffers from more severe seeing-induced calibration noise, and furthermore the 
target star WR 104 is also two magnitudes fainter here than at K (due to the rising IR 
excess). We therefore begin our discussion of the image morphology with images recovered 
from the CH4 filter in the K-band, and we later contrast the results from the other filters 
against this reference. 

The images of Figure [T] demonstrate the recovery of extended structures around the 
bright core with typical image dynamical ranges in excess of 100:1, although higher noise 
levels due to poor seeing did limit fidelity on occasion (e.g. JunOl and May04 epochs). The 
predominant feature at all epochs is a curved plume originating in a bright compact core, 
and wrapping around a full turn at which point there is a rapid fading in surface brightness 
to about 1% of the image peak. Beyond the bright first coil, at angular radii > 100 mas, 
there can often be seen further windings of a continuing spiral, or in some cases substantial 
segments of arcs, at a significance at or somewhat above (e.g. Apr98, Jun98) the noise level 
in each individual image. 

The overplotted dashed line on each image in Figure [T] gives the best-fit Archimedian 
spiral model, which is discussed in greater detail in the following section. Each image has 
been registered so that the mathematical center of the best-fit spiral (and not the brightest 
pixel) defines the coordinate origin. 



- 7- 



2.2. Archimedian spiral model 

An automated procedure was used to extract the coordinates of points lying along the 
crest of the bright spiral ridges from the images of Figure [D A model could then be fit to all 
11 epochs of data simultaneously. This model consisted of a uniformly-rotating Archimedian 
spiral, with free parameters: P the rotational period, W the angular expansion (wind) speed, 
% the inclination to the line-of-sight, the position angle of i, and 9q the position angle of 
the model at the first epoch (Apr98). There were two additional translational degrees of 
freedom associated with the registration of the model to the image at each epoch. 

With nearly 10 complete rotations of the structure covered by our data series, it was 
possible to obtain significantly more precise estimates for the model parameters than pre- 
viously published work. In particular, best-fit models constrained P to be 241.5±0.5 days, 
a value in very go od agreement with the previous best interferom etric value (243. 5±3 d; 



Tuthill et al.ll2002l ) and with the reported photometric period (241d; lKato et al.ll2002l ). The 



windspeed W was found to be 0. 28±0. 02 mas/day (or 102±7mas/yr), also in accord with 



prior findings (jTuthill et al.ll2002l ). Our best fit for i = 12°, although the 1-sigma uncer- 
tainty covers a range from 0° - 16°. Precise constraint of the inclination was hampered by 
the nearly face-on viewing angle to the system: 12° only results in a 2% distortion of the 
aspect ratio. Sky orientation parameters were found to be = 84° ±15° and 9q = 269° ±1°. 

Despite the fact that the model was only fit to ridge points extracted from the bright 
inner winding of the spiral, it can be seen from Figure [1] that it also provides a relatively 
good representation of the structure in the fainter outer regions. This implies that the second 
winding of the spiral follows closely the expected trajectory from simple expansion, and that 
the plume remains relatively well defined and compact as it evolves. 

A key finding of the modelling study was that the mathematical origin of the best-fit 
spiral was always systematically offset from the center of the bright core of flux in the images 
in a highly consistent fashion, as can be seen from Figured! This "standoff distance" between 
the spiral center and brightest pixel was found to be 13.3 ± 1.4 mas which amounted to a 
rotational displacement downstream along the locus of the spiral of 85° ± 14°. Implications 
for the physical model of the dusty colliding-wind binary from the modelling and imaging 
presented here are discussed further in later sections. 



2.3. H and L-band images 



A series of 11 images recovered from H-filter data are presented in Figure[2l Although 
image fidelity is somewhat lower than for the CH4 filter, the inner bright spiral structure is 
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well displayed at all epochs. Similarly, Figure [3] shows images recovered from the PAHCS 
filter from the short-wavelength side of L-band. Data were only collected in this band for 
the first six epochs of our study (see Tabled]), and the resolution obtained is lower than 
the other filters due to the longer wavelength. The spiral structure of the inner coil is again 
readily apparent. 

Archimedian spiral models were fitted to all H and PAHCS images following the proce- 
dure described above. Model parameters such as rotational period and winding angle were 
found to be in accord with the results from the CH4 filter, however the lower quality of 
the image data meant that the fits did not significantly advance our knowledge of these 
properties. 

A more profitable strategy for interpretation of the H and PAHCS images was to fix 
the spiral model parameters P, W, i and <fr to their values obtained above, and to examine 
the new data for evidence of systematic changes from this baseline model. One relatively 
subtle change found was in the mean rotational orientation of the spiral (model parameter 
6q). It was found that the mean rotational displacement of models fitting the H images were 
retarded by 6° (counterclockwise) while the PAHCS images were found to be advanced by 4° 
(clockwise) compared to the CH4 reference model. This effect was seen in two independent 
data-fitting approaches, the first used a similar strategy to that above in fitting to points 
along ridge crests, while the second entailed the generation of an average CH4 image which 
was then rotated and fit to each trial H/PAHCS image. The latter procedure was found to 
be more robust as it was less sensitive to resolution effects which might blur out the exact 
location of the ridge crest. 

A strong note of caution needs to be emphasized here because these relatively small 
changes in rotation are obtained from comparison between images with differing angular res- 
olutions. In order to investigate the effect of the decline in image resolution going to longer 
wavelengths H-CH4-PAHCS, simulations were performed in which images were artificially 
smoothed, and models fit to determine if there were systematic changes in the recovered 
parameters. Indeed, we found clear and systematic changes of best fit rotation angle, mea- 
sured to be about 5° in a counterclockwise sense when the simulated resolution was degraded 
to imitate the transitions H-CH4 and CH4-PAHCS. This makes intuitive sense: lowering 
the resolution tends to blur the center together, and the projecting arm of the spiral ap- 
pears to rotate as more of it is lost to the expanding core. However what is interesting is 
that this resolution-induced rotational bias is in the opposite sense to the observed clockwise 
rotations of 6° H-CH4 and 4° CH4-PAHCS. It is therefore likely that the true magnitude 
of this rotation with observing wavelength is underestimated here by a factor of two due 
to the bias counteracting the signal. However, the fitting of any models to images at the 
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resolution limit is challenging, and spirals perhaps even more so as due to their self-similar 
nature. The confirmation and unveiling of detailed substructures within the WR 104 plume 
must therefore await more powerful telescopes or imaging arrays with still higher angular 
resolution. 



2.4. Lightcurves 

Although our imaging experiment was not specifically designed to yield flux measure- 
ments, it was possible to extract photometry within our three near-infrared filters to within 
~ 0.1 — 0.2 mag, although for some nights there were light clouds which compromised the 
measurement accuracy. These measurements are presented in Figure HJ and show the infrared 
fluxes to be fairly constant in the H and CH4 filters (to within a few tenths of a magnitude). 
A larger variation is seen in the PAHCS filter, although it lacks any clear sinusoidal signature 
with phase and may arise from higher noise levels due to atmospheric contamination (this 
filter lies near the edge of the infrared L band window). 



3. Basic Geometry and Physical Properties 



In order to illustrate the physical properties of the WR 104 pinwheel system in the 
discussions that follow, a cartoon of the basic geometry is given in Figure [5j A pole-on view 
onto a 20° opening half-angle spiral shock is depicted, together with key elements of the 
system such as the central binary star, the optically thin inner region of the plume and the 
dusty first and second coils. 

If we are able to measure the physical speed of the dust plume, then this may be 
combined with the apparent angular velocity to yield the distance to the system. Following 
Tut hill et al.l (119991 ), if we identify the termina l WR windspeed obtained spectroscopically 
as Voo = 1 220km/s ( Howarth fc Schmut zj Il992l ) then we may combine this with our proper 
motion of 0.28±0. 02 mas/day to yield a distance to the system of 2.6±0.7kpc. By far the 
dominant error is from the ~25% un certainty in the windspeed fou nd by comparing the 
results of different line profile studies (IRochowicz fc Niedzielskil 119951 ) (velocities as high as 



1 600km/s have been reported for WR 104 ([Torres et al 



198 



3))- 



The distance derived in this fashion is somewhat larger tha n the 1.6 kpc estimate from 
possi ble membership of Sg r OBI (jLundstrom Sz Stenholm 1984 ); a discrepancy which moti- 
vated lHarries et al.l (120041 . hereafter HMSK) to question the validity of this method. These 
authors preferred the smaller distance scale on the basis that the WR wind may be subject 
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to radiative braking at the OB star (IGayley et al.lll997l ). and thus the terminal wind velocity 



may not be representative of the bulk motion of the observed dust plume. It is instructive 
to pursue this controversy further, for it illuminates some of the unique and interesting 
properties of these Pinwheel systems. 

Let us consider a very simple picture in which we distinguish only 3 populations of gas 
in the colliding wind system: the free WR wind, the free OB wind, and the material around 
the shock surface which will contain both, but with modified density and velocity. We will 
furthermore, for the present, restrict our attention to the hydrodynamic treatment of the 
inner shock with opening angle given by HMSK. This ignores possible modifications to the 
shock geometry due to radiative braking, and we devote the Section 0] to a discussion of 
these effects. 

For an ordinary (non-rotating) bowshock system, these populations of gas interact 
strongly in the heart of the system but at large distances are able to flow more-or-less 
unimpeded for they occupy distinct regions of space. As the dust is known to form around 
the colliding-wind interface, the arguments of HMSK would be correct in suggesting a lower 
velocity for the material in the dust plume from the various braking mechanisms (radiative, 
shocks). 

However, this picture breaks down for the Pinwheel systems when we introduce the 
rotation of the plume driven by the binary orbit. Recent three-dimensional hydrodynamical 
simulations of colliding-winds in binary systems which include the coriolis eff ect have shown 



strong departures from the simple non- rotating case (ILemaster et al.1 120071 ). Asymmetric 
shock strengths on the leading and trailing edges of the shock spiral were demonstrated. 
Although these simulations were mostly concerned with the inner regions of colliding wind 
systems, there will also be important consequences at large radii which need to be considered 
for the case of WR 104's dust plume. 

The resultant geometry will be much more complicated because the OB wind will be- 
come thoroughly wrapped and entangled within the spherically-expanding WR wind. The 
interpenetration of the two winds leads to a very different outcome than for the non-rotating 
case. It is no longer possible for the distinct populations of gas to travel at different speeds 
without interfering, and in the orbital plane containing the shock, at large radii all popula- 
tions of gas must collide and eventually come to a common velocity. The two key questions 
we need to answer to understand the form of the spiral plume are (1) what is this final 
velocity Vp, and (2) how long will it take until the material in the plume is moving at Vpl 

We may answer (1) with a quantitative discussion of the wind kinetic energy. Following 
HMSK, for the WR component we adopt a mass-loss rate M = 3 x lO~ 5 M /yr, and a 
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windspeed V^wr — 1220km/s, while for the OB star we have M = 6 x 10 8 M Q /yr and 
Voqob ~ 2 000 km/ s. Th e open ing; angle for the shock can now be calculated from the formula 



given m 



Eichler fc Usovl (119931 ) to be 9 = 17.4° half-angle (rounded to 20° half or 40° full 
opening angle in HMSK). Note that this explicitly ignores any widening of the opening angle 
due to radiative braking (see Section H]), but we proceed with this value for the present to 
ensure consistency with HMSK. At large radii, the solid-angle subtended by the shock cone 
will be Q = 27r(l — cosO) = 0.03 sr, so that in the non-rotating case, about 3% of the WR 
wind and 97% of the OB wind will participate in the collision. We define Rbowshock to be the 
volume ratio of the WR wind involved in the initial bowshock. 

In the rotating case, however, a far larger fraction of the WR wind will eventually 
collide due to the entanglement of the winds. Winds along the two polar axes, with solid 
angle Q p = 2n(l — cos(ir / 2 — 9)) = 4.1 sr from the WR can propagate unimpeded. However all 
gas within an equatorial angular band within ±9 from the plane of the orbit must eventually 
collide with the plume - for our geometry of 9 = 20° this corresponds to 34% of the total 
WR wind. Of this 34%, it is crucial to point out that only Rbowshock = 3% will undergo a 
'prompt' collision at or just downstream of the shock stagnation point, the remaining 31% 
will only collide at larger radii as the spiral windings gradually entangle the winds. 

If we assume, as seems reasonable for a relatively wide colliding-wind binary, that the 
shock is mostly adiabatic with little radiative energy loss. The kinetic energy flux of the 
asymptotic flow for the entangled winds can be found: 



-M F Vp = -RwrM W rV^ )WR + -RobMobV^qb (1) 

Here, Rwr and Rob are the fractional volumes of the WR and OB winds that are 
involved in the collision; from the geometrical arguments above we take Rwr — 0.34 and 
Rob = 1-0. Substituting values into this equation we calculate that the fast (but tenuous) 
OB wind will modify the final Vp in the equatorial plane from 1220km/s to 1226km/s: an 
insignificant change. 

The model may be further refined by adding a third population of radiatively braked 
material, as suggested by HMSK. Let us suppose (as an upper limit) that all WR wind 
impinging upon the 9 = 20° conical shock (3% of the total) is so braked. A preliminary 
estimate of 800km/s for Vp the streaming velocity of material in the shock cone has been 
obtained by fitting to the C III 5696 A line profile (G. Hill et al., 2007 in preparation). If 
we add this third population of gas to Equation [TJ, we obtain Vp =1171 km/s: a decrement 
of only 4% on Voo^/r- We therefore answer part (1) of our question above that Vp = V^^wr 
with only very minor adjustment, confirming that the momentum-dominant WR wind will 



- 12 - 



overwhelm other wind populations. 

The second part of the question, finding how long will it take for the WR wind to 
predominate, has a bearing on the shape of the inner spiral: it will be distorted away from 
Archimedian if the matter is accelerating. A rough estimate can be obtained from a simple 
"mass loading" argument applied in the equatorial plane of the system. Let us suppose that 
in the wake immediately behind the OB star, a mass Rbowshock x Mwr of matter in the plume 
has a velocity Vp <C V^wr- If material were to continue on this trajectory, after time t it 
would be overtaken by a mass of the faster wind equal to 

m = {l--^—)^R WR M WR (2) 

V oo,WR. r 

where P is the orbital period. We may solve this equation for t/P to give the fraction of an 
orbital period in which a mass of material has collided with the plume from the free WR 
wind which is equal to that already present from the initial bowshock. Using values from 
above, we find that 0.1 of an orbital period is sufficient to double the piled- up mass at the 
inner wall (see Figure EJ), and by the time matter has reached the dust formation zone at a 
quarter turn, it will be strongly (~70%) momentum-dominated by the V^wr wind. 

This analysis confirms the general picture that the spiral plume in this system is rapidly 
and comprehensively momentum- coupled to the strong WR wind. The implied distance 
should therefore be 2.6±0.7kpc at which the displacement between successive spiral coils is 
170 AU, and the excellent match to an Archimedian spiral confirmed. It is very interesting 
to speculate that the spiral geometry, which allows a far larger interaction between the winds 
and generates active shocks at much greater radii than a non-rotating system, may also help 
play a role in the physics of dust nucleation in these systems. Indeed, due to the differing 
velocities of gasses in the region immediately behind the wake, material will be piled up 
both at the inner and outer walls (see Figure Numerical simulation in three dimensions 
of this fascinating 'twist' to a conventional bowsh ock geometry, extending the coriolis force 



hydrodynamic treatment of iLemaster et al.l (120071 ) to large radii, may prove a valuable guide 



for further progress in understanding these structures. 



4. The Potential Role of Radiative Braking in WR 104 

For WR+OB binary systems in which the momentum of the Wolf-Rayet wind sub- 
stantially exceeds that of the OB star, the WR wind can penetrate very deeply into the 
acceleration region of the OB-star wind, and in such cases it becomes relevant to consider 
the role of the OB-star light in providing the momentum balance against the WR wind. 
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The analysis by Gayley, Owocki, & Cranmer (1997, hereafter GOC) suggests that this can 
lead to a "sudden radiative braking" of the WR wind that can hold the stagnation point 
along the line of centers farther off from the OB star. This could also modify the angle for 
the global wind interaction cone, with ramifications for the observations of dust formation. 
GOC further provide scaling relations to identify whether such radiative braking is likely 
to be of importance for any given system, and even provide a diagram (their figure 5) to 
classify several of the known close WR+O binary systems. Let us now apply these criterion 
to ascertain the potential role of radiative braking for WR 104. 

Note first that the systems identified by GOC as likely to have radiative braking are 
mostly relatively close binaries with periods on the order of days or weeks, and separations 
of only a few OB-star radii. For wider systems, spherical expansion dilutes the WR wind 
more strongly in the region around the OB-star, making it easier for the OB-star wind to 
maintain a standard wind- wind ram pressure balance. But if the WR/OB wind momentum 
ratio, 

PWR/0B = Mol^' <3) 
is sufficiently high, then even in a relatively wide system, such as WR 104, the standard 
ram pressure balance could potentially be supplanted by radiative braking. The model wind 
parameters for WR 104, from Section [31 do indeed imply an extreme momentum ratio, 
Pwr/ob — 305, and so even the relatively wide separation of D = 2.4 AU, or d = D / R Q b — 
50 OB-star radii, might still induce radiative braking prior to achieving the hydrodynamical 
stagnation of the winds. 



4.1. Estimating the braking radius in WR 104 

The following discussion makes use of a number of parameters that are described in 
greater detail in GOC; for convenience here we provide a brief glossary introducing significant 
quantities. Radiative acceleration is assumed proportional to r] which is found from the 
effective opacity due to the integra l of th e forest of spectral lines in the wind. The parameter 



a was introduced by ICastor et al.l (119751 ) to describe how sensitively the line force increases 
with wind acceleration. GOC found it convenient to define normalized quantities for radiative 
braking separation d r \, and momentum ratio P r b which could be used to scale the actual 
separation and momentum ratio to gauge the significance of braking in individual cases (see 
GOC for detailed derivations). 

To locate WR 104 on the GOC figure 5 classification diagram, this wind momentum 
ratio Pwr/ob an d separation d must be further scaled by associated radiative braking values, 
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P r b and d r b- To determine those, we first need to evaluate the parameter 77, given by GOC 
eqn. (29), which unfortunately includes an undefined parameter v* (for that parameter, see 
Gayley, Owocki, & Cranmer 1996). Reconstruction of the GOC analysis leads to 

1 ( Lob \ ^ 1 °^ 2GM WR 4 / L OB \ 2 ( V eS c,WR \ 2 R WR 

V 1 : a«/0-<*){l + a) \L WR ) V^ R R OB 1 1 3 \L WR ) \ V WR J Rob' U 

where the latter equality applies to the simple case that the power index a = 1/2, and also 
recasts the scaling in terms of the ratio of WR wind speed to surface escape speed. Note 
that, in addition to replacing the undefined t>*, this formulation has the advantage of being 
manifestly independent of the separation D, as required by GOC eqns. (11), (12), and (17). 

Applying the wind parameters from Section [3], we find (assuming a = 1/2) that for 
WR 104, 77 = 1.65. Plugging this into GOC eqn. (30) then leads to d rb = 2.1. For a 
standard velocity index (3 — 1, we then further find P r b = 0.642. Using these with the above 
values for Pwr/ob an d d, we then find a scaled momentum ratio 

P = = — = 475, (5) 

P rb 0.642 ' v ; 



and scaled separation 



? d 50 n * 

d = — = — = 24. (6) 
d r b 2 . 1 



In the logP vs. logd plane, systems with radiative braking are defined by the triangular 
region with both log d > and log P > 2 log d. The former ensures that the momentum 
of the OB-star light is sufficient to keep the WR wind from impacting the OB-star surface, 
and is easily satisfied in this case. The latter condition implies that there can be no normal 
ram pressure balance between the winds, and thus that radiative braking must stop the WR 
wind; since for the above parameters, we have P = 475 > 407 = d 2 , implying that this is 
condition is also satisfied, although more marginally. If, for example, the estimated wind 
momentum ratio were reduced by more than about 15%, or the separation were increased 
by more than about 7%, then the system would no longer satisfy this latter condition. In 
those normal ram-pressure balance is still allowed, but even then radiative braking 

can play a role if the braking radius is outside the wind-wind balance radius. 

From GOC eqn. (18), we find that the braking radius for the above standard parameter 
set is 

r b = x b D = 1 + ( ^ )1/3 = 0.254 D = 12.1 Rob. (7) 

For comparison, for the simple case of constant flow speed (i.e. velocity power index (3 = 0), 
ram pressure balance in the absence of direct radiative forces occurs at a radius (cf. GOC 
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eqn. (2)) 

r s = x s D = - D = 0.054 D = 2.6 R OB ■ (8) 

1 + V "WR/OB 

For such a purely hydrodynamical balance, a global momentum analysis (Canto, Raga, 
and Wilkin 1996; Gayley & Owocki in preparation) yields a transcendental relation between 
the cone-opening half-angle 9 S and the momentum ratio Pwr/ob, 

tan0 s -0 s = -. (9) 

^WR/OB — J- 

For large momentum ratio Pwr/ob 3> 1, use of a small-angle approximation for tan# s = 
9 S + 9^/3 gives the asymptotic form, 

/ 3tt \ 1/3 

S = ; Pwr/ob > I- (10) 

V rWR/OB J 

To recover 9 S = 90° for the symmetric momentum case Pwr/ob = 1, we can write a more 
generally applicable approximate solution, accurate to within about 6% for all Pwr/o, 

Q s = 1/3 (11) 

31/90 + P WR iq B 

where using 121 = (180/7r) (Stt) 1 ^ 3 gives this angle in degrees. 

For our estimate of Pwr/ob = 305 for WR 104, we then obtain 9 S = 18° , in good 
agreem ent with the value from Section [3] derived from the heuristic formula of lEichler fc Usov 
(jl993h . 



4.2. Radiative force influence on the shock-cone angle 

An important question is then, by how much can the effects of radiative forces from 
the OB star potentially widen this shock cone? Only a complete simulation can answer this 
definitively, but we can make conceptual progress using heuristic approximations to attempt 
to better constrain the circumstances under which we should expect an important radiative 
modification of the bow shock geometry. Two approaches that borrow from the existing 
sudden radiative braking analysis can be termed the "Global Momentum Augmentation" 
(GMA) and the "Radiative Wall" (RW) approach. Each gives similar results for the opening 
angle, and so here we briefly outline just the former, deferring to Gayley & Owocki (in 
preparation) for a detailed discussion of both approaches. 
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In the GMA approach, the radiation interaction with the incident WR wind is assumed 
to augment the OB-star wind momentum flux by a globally constant factor, set by the 
requiring that the stagnation point of the combined effective momentum fluxes equal the 
braking radius in eqn. (j7j). We write this factor as jt, where 7 is the ratio of the OB-star 
radiative momentum flux to its wind momentum flux, Lob/ MobVobc, and r is an optical- 
depth-like parameter that simulates the coupling of the OB radiative momentum flux to the 
WR wind. In analogy to the global analysis cited above (Gayley & Owocki, in preparation), 
we find for the large momentum ratio limit, 

0. = 121 f^^y 73 • (12) 

In this scenario radiative braking would alter the shock geometry whenever 7T is ap- 
preciable, and since 7 in WR 104 is expected to be large, this is not a terribly restrictive 
requirement for the coupling parameter r. Since the braking radius scales roughly with 
rj^ 1 / 3 when braking is important, the condition that the effective momentum ratio used in 
r s be reduced enough to achieve r r b = r s (see eqs. [7] and [8]) is tantamount to requiring 

^ / vRob\ 2/3 p 

IT = I — — J PwR/OB • (13) 

This implies a shock opening angle 

e, 2 121 (^f) 2 " , (14) 

which for 77 = 1.65 and D/Rob = 46.5 evaluates to 6 S = 58°, representing a factor ~ 3.2 
increase over the purely hydrodynamic value 9 S = 18°. This may be a significant overesti- 
mate, because of the optimistic assumptions made about the radiative momentum coupling, 
but it demonstrates that radiative braking could indeed increase the shock cone opening if 
the WR wind opacity continues to be high throughout the flow. 

A similar result can be reached through the radiative wall approach. Although both 
approximations are too heuristic to be considered quantitatively reliable, they do suggest 
that strong radiative braking for the assumed parameters in WR 104 should be expected to 
widen the bow shock angle considerably, possibly by a factor of 2 or more when compared 
to the purely hydrodynamic result 9 S = 18°, as long as the incident WR wind maintains 
the same high levels of opacity for the OB starlight that it must have had to have been 
radiatively driven from the WR star. This assumption is currently of unknown validity, and 
indeed if we take the opposite limit, and say that the line-opacity is just what is needed 
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to drive the relatively weak OB-star wind, then radiative braking will be much weaker and 
would not likely alter the shock-cone geometry in any measurable way. 

To address intermediate values of the effective WR wind opacity, recalling that 77, for 
fixed a, is simply proportional to the integrated spectrum- weighted line opacity distribution 
with no regard to optical depth (i.e., self- shadowing) corrections. Thus one may either make 
assumptions about this line opacity to predict rj, or one may use observations of the bow 
shock cone to infer r\ and reason to conclusions about the line distribution. Further analysis 
will be required to constrain these expectations more quantitatively, but it seems at least 
plausible that radiative braking may be important in the wind-wind interaction and shock- 
cone geometry of WR 104, if the line opacity of the WR wind interacting with the OB 
starlight is substantially enhanced above what is needed to drive the OB-star wind itself. 
Reversing the logic, it may be concluded that evidence for a significant widening of the 
bow-shock cone relative to what may be explained hydrodynamically may be interpreted as 
evidence that the WR wind does indeed carry with it an intrinsically high opacity when it 
impinges into the region of strong OB-star radiative flux. 

Our images of the WR 104 plume are able to rule out some of the more extreme enlarge- 
ments to the cone opening angle in the above scenarios. Fits to the data by HMSK found 
that model cone half-angles of 40° or greater (twice the hydrostatic value) were rejected. 
With still higher angular resolutions than those obtained here, it should be possible to get 
a direct measurement of the cone opening angle from the images, revealing the role played 
by radiative braking on the shock geometry. 

5. The Profile of the Spiral Plume 

Using the excellent fits provided by the theoretical spiral locus to the observed plume of 
emission from Figures [H El & El an investigation of the brightness profile along the length of 
the dust plume was performed. The flux, summed over a band following the model spiral, was 
accumulated over all images recovered. Average profiles, plotted as a function of angular 
displacement along the spiral, are given in Figure [6] for all three filter bandpasses. Also 
overplotted is the same profile along the spiral calculated from the synthetic radiative-transfer 
images of HMSK. All experimental curves follow a similar form (with some departures for 
PAHCS which we attribute to the systematically lower resolution) exhibiting a steep rise 
in flux from the spiral center to a point about 90° downstream followed by a rapid drop as 
the bright core of the system is traversed. There now follows a region of roughly monotonic 
decline between about 150 — 400°, at which point there is a very sharp drop in brightness 
flattening to level of less than 1% of the peak for the second coil (CH4 filter only). 
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The initial sharp rise in flux seen in Figure [6] gives another way to visualize the "standoff 
distance" between the spiral center and brightest pixel, measured earlier to be 13.3 ±1.4 mas 
or 85° ± 14° downstream along the spiral. Given the 241.5 d period, this corresponds to a 
kinematic age for the material of 57 ± 10 d. This value is in remarkable agreement with the 
value of ~ 58-87 d for the delay between periastron passage and the peak in t he dust produc- 



tion in the well-studied episo dic colliding-wind binary dust producer WR 140 (IWilliams et al. 



19901 ; iMarchenko et al.l 120031 ). The finding that the peak in dust production occurs at some 
distance downstream from the colliding-wind stagnation point (~30 AU for WR 104) points 
to a more complicated picture for dust for mation phys ics in colliding-wind binaries than the 



simple shock compression/cooling mo dels (lUsoy Il991l) . A number of avenues for additional 



physics have been presented (see also ICrowther 



200 



J, with one likely element being escape 



from the dust-hostile immediate circumstellar radiation environment. Chemical seeding of 
the WR wind (predominantly helium) with hydr ogen from the O B star at the contact surface 
may help to enhance dust formation pathways (ILe Teuffl 120021 ). although such mixing may 
not be very rapid. Shielding from the intense radiation field does appear to be a key ingredi- 
ent, and sets the playing field for sharply nonlinear grain growth in which small clumps may 
self-shadow larger areas downstream in a runaway process. Recent numerical simulations of 
colliding flows have shown that matter is squeezed into thin high- density shells threaded by 
still higher density filaments as a result of supersonic turbulence ( IWalder &: Folinil |2000| ) . It 
was suggested these conditions can persist for a considerable time in the flow before dissipa- 
tion by internal shocks and vortex cascades. Further to these earlier ideas, we now also add 
the idea of wind entanglement which can generate significant shock activity well downstream 
of the nose, as discussed in earlier sections. Whatever mix of these or other elements may 
go into the next generation of colliding-wind dust formation models, there certainly seems 
scope to address the outstanding issue of whether the dusty- WC phenomena has a necessar y 
link to binarity (IMonnier et al.lll999l ; IWilliams fc van der Huchtll2000l ; iMonnier et al.l 120071 ). 
or whether novel du st formation pathways in an isolate d wind may also be important in 
some circumstances (IZubkolll998l ; ICherchneff et al.ll2000l ). 

After the peak at 85°, the profiles of Figure [6] exhibit a steady decline (but note the 
logarithmic y-axis). This is in accord with the relatively flat temperature profile found 
in models of the centrally-illuminated dust spiral models of HMSK, although the bright 
peak near the heart of the system is not reproduced in the simulations. A fairly detailed 
comparison between the radiative transfer results and published masking images of WR 104 
can be found in HMSK, although our new results here highlight some discrepancies which 
we hope will motivate future refinements. The flux level of the second winding of the spiral 
is well reproduced at the beginning (480°-600°) but the secondary bright knot in the model 
tail seen at 2 complete turns was not observed in the data. Although a detailed comparison 
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is beyond the scope of the present paper, one possible implication is that the optical depth 
of the dust in the model is not high enough in the core of the system. Increasing it may help 
to generate the bright peak at 85°, and at the same time decrease the illumination on the 
second coil therefore suppressing the spurious feature at 720°. 

The very sharp decline in flux between 400-480° is a result of the shadow cast by 
the first winding of the spiral onto the second, causing a rapid drop in temperature for the 
eclipsed material. T his feature has been commented upon since the discovery of the pinwheels 



( ITuthill et al.lll999l ). and was also reproduced in radiative transfer models (HMSK). What is 
particularly interesting is that the terminator is found at an angular ordinate of ~ 440° along 
the spiral, or about 80° beyond one complete revolution. This extra angular displacement 
(easily visualized in the figure as the full revolution 360° location has been marked) gives 
further independent confirmation for the existence and size of the optically thin inner cavity 
surrounding the central stars. We refer the reader back to Figure [5] for an easily understood 
and intuitive depiction of the angular locations of the inner optically thin plume, and the 
geometry of the shadow cast upon the second coil by the first. 

Also of interest in helping to constrain future generations of dust plume models are the 
multi- wavelength observations of the spiral given in Figures [U|2]& [31 In particular, there was 
a systematic rotational displacement found with observing wavelength with a (lower limit) 
position angle of 6° from H-CH4 and 4° from CH4-PAHCS (true values may be a factor 
of ~2 larger - see discussion in Section 12.31) . This implies that the transverse temperature 
profile of the plume can be explored by our high resolution techniques. Short wavelength 
observations will preferentially weight the hottest material on the inner edge of the plume 
facing the central stars, while longer wavelengths will give a better indicator of the bulk 
concentration of warm dust. Again, we must defer detailed modelling of these effects to a 
future paper, but note that preferential radiative heating of the inner wall was predicted in 
models of HMSK. An alternate scenario is that of asymmetric heating by shocks of differing 
strength at th e inner and outer walls, as recently demonstrated in the rotating hydrodynamic 



simulations of iLemaster et all (120071 ). 



6. Geometrical Implications of the Plume Structure 



A key finding from earlier st udies of the WR 104 pinw heel has been the l ow eccentricity 



of th e orbit of the central stars (iTuthill et al.lll999l . 120021 ). As noted earlier (jMonnier et al. 



19991 ). this finding may be a "smoking gun" pointing to previous episodes of tidal circu- 



larization and likely Roche- lobe interactions with possible mass transfer and/or envelope 
stripping. Such an event provides one of the two possible pathways to the creation of a 
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Wolf-Ravet ( 


Paczvriski 


1967 




(Chiosi et al. 


|l978 


). With our 



19781 ). With our long time-sequence of observations sampling all phases of the 
orbit, we are able to put tight constraints on the allowed levels of orbital eccentricity in the 
case of WR 104. Although it would be possible to fit eccentric spiral models directly to the 
data, we favored a more indirect approach as for small to moderate levels of eccentricity, the 
departure in shape from an Archimedian spiral was subtle. The most easily detectable effect 
that an eccentric orbit has in perturbing the rotating spiral model away from its ideal form 
will be in departures from a constant rate of rotation. 

We have therefore performed a careful fit of the best-fit orientation of the spiral structure 
at each epoch, and compared the values obtained with the model of a simple, constant angular 
velocity. This was done in two different ways which gave similar results: firstly by fitting 
to points extracted from the crest of the bright spiral ridge, and secondly by fitting each 
full image to a mean template. The uniformity of the angular rotation rate extracted from 
the real data, which is given in in Figure [7J was then compared with results from model 
simulations. A family of synthetically-generated spiral images of varying eccentricity were 
fit in an identical fashion to the real data, allowing us to derive quantitative limits on the 
levels of eccentricity present in the WR 104 orbit. Figure [7] shows the rotational uniformity 
from the data and examples from two of the synthetic eccentric model images. The rotation 
is seen to be highly uniform, ruling out all but very modest eccentricities. A 2<j upper limit 
of e < 0.06 was obtained in this fashion. 

Having thus established the properties of the best-fit uniform Archimedian spiral model, 
it was then possible to generate a stacked composite image averaging over all epochs of data 
taken in the CH4 filter. Individual data frames were corrected for inclination on the sky 
(simulating a perfectly face-on view) and de-rotated to coincide with the first (Apr98) epoch 
before being co-added. A composite image, displaying significantly greater dynamic range 
than any of the single-epoch images from Figure [TJ is given in Figure [83 Variations in 
morphology between different epochs beyond those expected from the rotating spiral model 
are minor, and all appear to be within the noise level of the image recovery process. We 
therefore hope that this image, as a mean over all epochs, may provide a resource to modellers 
and others interested in the best possible snapshot of the geometry of these pinwheel systems. 

The flat infrared lightcurve s of Figure H] d iscussed above affirm the status of WR 104 



as a "constant" dust-producer (j Williams! 1 199 71 ). yet stand in contrast to the large 2.7mag 



"q uasi-periodic" vari ability (periodic component > 1 mag) in the visible lightcurves reported 



by iKato et al.l (120021 ) . With a narrow cone half-angle (20°; HMSK), face-on orientation to 
the line of sight and very low levels of orbital eccentricity, the spiral dust plume gives no 
obvious mechanism for modulating the visible flux from this system. The lack of detectable 
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infrared variation argues for steady and constant dust production at all orbital phases, with 
minimal effects from any inclination to the the line of sight. Visible light scattered from the 
dust plume may be more strongly modulated by rotation and inclination effects, however 
it is difficult to see how variations in the V lightcurve as large as > 1 mag could arise from 
this. Possible scenarios include opacity effects from gas (or possibly small quantities of 
dust) associated with the innermost regions of the shock cone which modulate the light 
of the OB star, line-of-sight (polar-axis) dust creation independent of the colliding wind 
mechanism, or some geometrical effect from a sh arp opacity gradi ent possibly created at an 
earlier phase (further discussion can be found in iKato et alJEooS ). 



7. WR Systems as Supernova Progenitors 

As a final point of discussion, we also explore the implications of future evolution of this 
system with particular regard to the face-on viewing angle to Earth. The formal Archmiedian 
spiral fits yielded an inclination i = 12°, but for such a small angle to the line-of-sight the 
uncertainty was large and fully consistent with the range from 0° - 16°. A further source of 
uncertainty arises from the limitations of the simple two-dimensional spiral model used. In 
reality, the pinwheel is of course a three-dimensional structure which may present us with 
more complicated opacity effects such as limb-brightening. If we consider a single circular 
cross-section through the conical shock, as this is carried downstream by the WR wind, it will 
become inflated and distorted by the spherical expansion into an elongated "D" . Eventually 
at large radii, the plume shape will resemble an equatorial band wrapped around a sphere, 
subtending an angle equal to that of the shock cone. This underlying spherical geometry 
gives rise to a fundamental insensitivity to the line-of-sight inclination: modest changes of 
tilt for such a structure result in only a very small or zero change in observed aspect ratio. 
For these reasons, our upper limit of 16° may be an underestimate. Firmer constraint of the 
inclination from the data presented here would require detailed radiative transfer modelling 
with varying lines of sight onto a three-dimensional structure. 

Despite these quantitative caveats, it remains apparent that in WR 104, Earth has an 
approximately pole-on view onto a relatively nearby, massive WR binary system. Given 
the prehistory of mass transfer and envelope-stripping suggested above, it therefore seems 
likely that Earth also lies in the polar direction with respect to the spin axis of both WR 



Zahn 


1977; 


Tassoul 


1987) 



core-collapse supernovae, it is interesting to speculate on how our "privileged" polar vantage 
point may affect our experience of this event, which should happen within a timescale of 
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sever al hundred thousand years for the WR component (IMaederl Il98ll ; IChiosi fc Maeder 
19861 ). It is becoming more accepted that all supernova explosions probably exhibit some 
preferred axis aligned with the progenitor stellar spin. As one of the closest imminent WC- 
class Wolf-Rayets, an isotropic supernova in WR 104 at kiloparsec scale distances would 
no doubt put on an impressive sh ow, but any impact on Earth's biosphere is likely to be 
negligible (jEllis fc Schramm! 119951 ). However for a highly anisotropic explosion, the most 
extreme example of which would be a Ga mma Ray Burst, effects could be significant if 
Earth lies within the £ 12° opening angle ( Frail et al.l 2001 ) of the burst, even at 2kpc 
distances flMelott et aDbooi Ehomas et alil2005al lbl). 



Could WR 104 produce a GRB? Certainly a case could be made that the OB component 
of the binary system, which was the recipient of mass and angular-momentum transfer in the 
prior Roche-lobe overflow event, may meet the requirements for envelope mass an d angular 
momentum for many current models of GRB formation (see iPetrovic et al.ll2005l . for a dis- 
cussion of WR binaries as GRB progenitors). The Wolf-Rayet component, with WC spectral 
type, has arguably the shortest fuse, being in the last known stable phase before supernova. 
Although prior stripping may have lowered the envelope mass into a regime less favorable for 
conventional collapsar GRB's, there ar e also more exotic magnetar (millisec pul sar driven) 
GRB models which may prove viable (IGaensler et al.ll2005l ; iMetzger et al.l 120071 ). In short, 
this is a very active and rapidly-evolving field, and the uncertainties outnumber the firm con- 
clusions in arriving at a mass/metallicity/evolutionary portrait of likely GRB progenitors. 
Whether or not the WR 104 system will play host to a future GRB, it does seem clear that 
two pole-oriented supernova e will occur, and that energy and matter will be preferentially 
ejected along this axis (see iFryer fc Warren! 120041 . for a discussion of the effects of stellar 
rotation on SN explosions). Further observations to constrain the inclination of the orbital 
plane, which should be possible with spectroscopic monitoring, are therefore encouraged. 
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Fig. 1. — Images recovered from data in the CH4 filter over II observing epochs. Contour 
levels are .4, 1, 2, 5, 10, 50 % of the peak. North (Dec.) is up and east (R.A.) to the left with 
the image scale labelled in millseconds of arc. Images have been centered on the mathematical 
origin of the best-fit Archimedian spiral model (dashed line: see text for details). 
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Fig. 2. — Images recovered from data in the H filter over 11 observing epochs. See caption 
to Figure [1] for details. 




Fig. 3. — Images recovered from data in the PAHCS filter over 6 observing epochs. Contour 
levels are 1, 2, 5, 10, 20% of the peak. Otherwise similar to Figure [U 
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Fig. 4. — Photometry data for WR 104 obtained in our three filter bandpasses (key inset), 
folded with the 241.5 d period. 
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Fig. 5. — Schematic diagram of the key elements in the WR 104 system. The geometry 
depicted is that of a face-on view onto a 20° cone half-angle pinwheel. The axis of the 
system, and the line of shadow cast at 84° to the spiral origin (6 = 0), are indicated as 
dot-dashed lines. Key elements are labeled, and some, such as the binary separation, are 
not to scale (depicted here about a factor of ten too large). 
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Fig. 6. — Average brightness profile of the plume as a function of angular displacement along 
the spiral model locus. Data for our three filter bandpasses are given (linetypes indicated in 
th e key) together with a profile extracted from the numerical radiative transfer simulations 
of lHarries et al.l (120041 1 specific to our CH4 filter in the K band. Data for H and PAHCS are 
terminated shortly after the first full turn, at which point the signal-to-noise in the images 
is insufficient to discern real structure from noise. For convenience, intervals of 1 and 2 
complete revolutions are marked. 
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Fig. 7. — (Upper panel) The rotation of the pinwheel at each observing epoch from fits 
to data in Figure (TJ plotted as a function of the orbital phase. The dashed line gives a 
simple model of uniform angular velocity. (Lower panel) the residual between the data and 
the uniform model. Both panels have overplotted dotted lines giving examples of expected 
curves for elliptical orbits with eccentricity e = 0.1 and 0.2. 
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Fig. 8. — Stacked composite image of all epochs of data in the CH4 filter. Contour levels 
are .1, .2, .5, 1, 2, 5, 10, 20, 30, 70% of the peak. Individual images have been projected 
to a face-on viewing angle and derotated to the first Apr98 epoch according to the best-fit 
Archimedian spiral model prior to coadding. 
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Table 1. Journal of Interferometric Observations. 



Date 


JD-2450000 


Filter 


1998 Apr 14 


918 


H,CH4,PAHCS 


1998 Apr 15 


919 


T T /" 1 T T A T 1 ATT f~1 

H,CH4,PAHCS 


1998 Jun 04 


969 


H,CH4,PAHCS 


1998 Sep 29 


1086 


H,CH4,PAHCS 


1999 Apr 25 


1294 


H,CH4,PAHCS 


1999 Jul 29 


1389 


H,CH4,PAHCS 


2000 Jun 23 


1719 


H,CH4,PAHCS 


2001 Jun 11 


2072 


H,CH4 


2002 Jul 23 


2479 


H,CH4 


2003 May 12 


2772 


H,CH4 


2004 May 28 


3154 


H,CH4 


2004 Sep 03 


3154 


H,CH4 



Table 2. Properties of Filters 



Name Center Bandwidth 
(/im) (/xm) 



H 

CH4 
PAHCS 



1.656 
2.269 
3.083 



0.333 
0.155 
0.101 



